High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

常微分方程技术及其在固体力学计算中的应用

常微分方程(ODE——Ordinary Differential Equation)边值问题的最新计算求解技术的迅速发展推出了一批高质高效的通用软件,而工程中大量的ODE问题并非呈现为这些求解器(Solver)所接受的标准形式。然而,运用一些简单的ODE变换技巧可以将大量的不同类型的特殊问题转化为标准形式。本文列举了若干常用的变换技巧,并广泛地应用于各种固体力学问题的计算中,使大量的ODE问题在形式上得到统一,得以用标准的ODE Solver方便有效地求解。     

DIFFERENTIAL-ALGEBRAIC APPROACH TO COUPLED PROBLEMS OF DYNAMIC THERMOELASTICITY

An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction procedure of the original system of PDEs describing coupled thermomechanical behavior to a system of Differential Algebraic Equations (DAEs) where the stress-strain relationships are treated as algebraic equations. The resulting system of DAEs was then solved with a Backward Differentiation Formula (BDF) using a fully implicit algorithm. The described procedure was explained in detail, and its effectiveness was demonstrated on the solution of a transient uncoupled thermoelastic problem, for which an analytical solution is known, as well as on a fully coupled problem in the two-dimensional case.     

A numerical investigation of transient detonation in granulated material

A two-phase model based upon principles of continuum mixture theory is numerically solved to predict the evolution of detonation in a granulated reactive material. Shock to detonation transition (SDT) is considered whereby combustion is initiated due to compression of the material by a moving piston. In particular, this study demonstrates the existence of a SDT event which gives rise to a steady two-phase Chapman-Jouguet (CJ) detonation structure consisting of a single lead shock in the gas and an unshocked solid; this structure has previously been independently predicted by a steady-state theory. The unsteady model equations, which constitute a non-strictly hyperbolic system, are numerically solved using a modern high-resolution method. The numerical method is based on Godunov's method, and utilizes an approximate solution for the two-phase Riemann problem. Comparisons are made between numerical predictions and known theoretical results for 1) an inert two-phase shock tube problem, 2) an inert compaction wave structure, and 3) a reactive two-phase detonation structure; in all cases, good agreement exists. Received 4 August 1995 / Accepted 17 February 1996     

三阶非线性发展方程的线性化有理逼近方法

研究了三阶非线性发展方程的初边值问题的解。采用基于Sinc函数的微分求积法发展了线性化有理逼近方法。通常的配点法不适用于上述三阶问题的求解。本文把提出的方法用于求解KdV方程,取得了良好的效果。     

A relaxation Riemann solver for compressible two-phase flow with phase transition and surface tension

The dynamics of two-phase flows depend crucially on interfacial effects like surface tension and phase transition. A numerical method for compressible inviscid flows is proposed that accounts in particular for these two effects. The approach relies on the solution of Riemann-like problems across the interface that separates the liquid and the vapour phase. Since the analytical solutions of the Riemann problems are only known in particular cases an approximative Riemann solver for arbitrary settings is constructed. The approximative solutions rely on the relaxation technique.The local well-posedness of the approximative solver is proven. Finally we present numerical experiments for radially symmetric configurations that underline the reliability and efficiency of the numerical scheme.